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Abstract. We revisit the XMM-Newton observations of the dwarf nova OY Car taken in July 2000 which occured 
shortly after an outburst. Ramsay et al (2001a) found a prominent energy dependent modulation at a period of 
2240 sec: this modulation was only seen for ~l/3 of the observation duration. In our new analysis, we examine 
this time interval in greater detail. In addition to the 2240 sec period we find evidence for other periods, the most 
prominent being near 3500 sec. Both these modulations are most likely due to changes in photoelectric absorption 
over this period: this is supported by phase-resolved spectroscopy. This may indicate the presence of matter above 
the accretion disc or a presence of a magnetic accretion curtain. In this case the 2240 sec period could represent 
a spin period of the white dwarf and the 3500 sec period a beat period between the spin and orbital periods. We 
also model the B band and UV eclipse profiles and light curves using a new technique to map the spatial extent 
of the accretion disc. As a result we find that whilst the optical emission is dominated by both the emission close 
to the accretion disc boundary layer and the hot spot where the accretion stream hits the disc, the UV emission 
is mainly dominated by the inner disc/boundary layer only. 
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1. Introduction 

Dwarf novae are interacting stellar binary systems in 
which material gets transfered from a red dwarf secondary 
star onto a white dwarf via Roche lobe overflow. In the 
absence of a strong magnetic field this material forms an 
accretion disc around the white dwarf. Dwarf novae show 
outbursts which occur on week to month timescales. For 
systems which are viewed close to the plane of the binary 
system, eclipses of the white dwarf by the secondary star 
are seen every orbital period. These systems are important 
since they allow a detailed study of the accretion disc as 
the secondary gradually covers then uncovers the disc. 

One such eclipsing dwarf nova is OY Car. It has a bi- 
nary orbital period of 1.51 hrs and has been extensively 
studied at various wavelengths. Recently, it has been ob- 
served at X-ray and optical/UV energies using XMM- 
Newton. Ramsay et al (2001a) found that there was a 
quasi-stable modulation of the X-rays at ^2240 sec, which 
was most prominent at lower energies. They speculated 
that this was the due to disc material orbiting above the 
binary orbital plane or due to the spin period of the white 
dwarf. In the latter scenario, this would represent the first 
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time that the spin period of the white dwarf in a dwarf 
nova had been detected. 

The eclipse profile shows a sharp drop in X-rays and 
in the UV band. In X-rays the flux remains at a very low, 
but significant level after this rapid (20-30 sec) drop. In 
the UV and B bands, the flux decreases at a more gradual 
rate after the rapid drop since the accretion disc is still 
visible after the eclipse of the white dwarf. Wheatley & 
West (2003) further analysed the XMM-Newton data and 
found that the X-ray eclipse lasted 30±3sec and concluded 
the X-ray emission originated at the poles of the white 
dwarf. 

In contrast, in this paper we investigate the nature 
of the 2240 sec modulation in more detail by revisiting 
the XMM-Newton data. In particular we have carried out 
phased resolved spectroscopy over this modulation period. 
Further, we apply a new eclipse mapping technique to 
model the spatial distribution of the accretion disc using 
the OM XMM-Newton data. 

2. Observations 

XMM-Newton was launched in Dec 1999 and has the 
largest effective area of any imaging X-ray telescope. It 
has 3 medium spectral resolution CCD type cameras on- 
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Fig. 1. The 0.15-1.0 keV power spectrum with 95%, 99% 
and 99.9% confidence limits (three dashed lines) from our 
simulations (top) and a power spectrum produced using 
simulated fake data that contains the two strongest signals 
(2193 sec and 3510 sec) seen in the power spectrum. The 
shapes of 2193 sec and 3510 sec pulses are estimated using 
a simultaneous 2nd order fourier fit with two periods. 

board: two EPIC MOS detectors and one EPIC pn detec- 
tor. Since the EPIC pn detector has the higher sensitivity 
we only discuss data taken using this detector (Striider et 
al 2001). In addition, it has an 0.3m optical/UV imaging 
telescope which enables simultaneous X-ray optical/UV 
observations (Mason et al 2001). 

OY Car was observed during the performance verifi- 
cation phase of the mission on June 29-30th 2000. The 
details of the observation are reported in Ramsay et al 
(2001a, b). However, the main points to note are that these 
observations were made some 4 days after an optical out- 
burst. The system had already faded by V ^4 mag from its 
peak, but it was still relatively bright in X-rays (~1.1 ct/s 
in the EPIC pn). The exposure in the EPIC pn detector 
was 48ksec, it was configured in full frame mode and the 
medium filter was used. The optical data were obtained in 
the UVW1 filter (effective wavelength 2910 A). The data 
were processed using the XMM-Newton Science Analysis 
System v5.3.3. Additional optical data in B band was ob- 
tained in August 2000 when it was at a similar brightness 
to the June observation (Ramsay et al 2001a). 

3. X-ray Light Curve Analysis 

Ramsay et al (2001a) found evidence for a quasi-stable 
X-ray modulation on a period of ^2240 sec, which was 
strongest in soft X-rays. This modulation was observed 



Fig. 2. The power spectra at 0.15-1.0 keV (top), 1.0-4.0 
keV (middle) and 4.0-10.0 keV (bottom). 



for around 1/3 of the duration of the observation and it 
was speculated that this period was related to the spin 
period of the white dwarf. 

We have re-examined the light curves of OY Car: in 
particular we perform a timing analysis of that portion 
of the light curve in which the 2240 sec modulation is 
strongest (HJD 2451725.6548-2451725.8622) in 4 energy 
bands (0.15-10.0, 0.15-1.0, 1-4 and 4-10keV). (We re- 
moved the time of the eclipse from the light curves). We 
show the amplitude spectrum of the 0.15-1.0keV light 
curve in Figure ^ (top panel) . The most prominent peak 
has a period of 2193 sec: within the errors, this period 
is consistent with the 2240 sec period found by Ramsay 
et al (2001a). However, there is a second prominent peak 
at 3510 sec: this peak is still present if we pre-whiten the 
data on the 2193 sec period. We performed the same anal- 
ysis in the 2 more restricted energy bands (Figure |2J): we 
find that both the 2193 sec and 3510 sec periods are more 
prominent at softer energies. A peak was seen at 3510 sec 
in the amplitude spectrum shown in Ramsay et al (2001a), 
however, it is more prominent in this new study since we 
have restricted the analysis to the time interval in which 
the quasi-stable modulation is seen. In order to demon- 
strate the prominence of both of these periods, we show 
the data folded on either of the two periods, after they 
have been prewhitened using the other period (Figure 13). 

To determine the significance of these two periods we 
randomly populated the time bins with the observed count 
rates. We then generated 10000 power density spectra: 
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Fig. 3. The 0.15-1.0 keV light curve prewhitened on the 
3510 sec period and then folded on the 2193 sec period 
(top panel), and vice- versa (bottom panel). 



the 95%, 99% and 99.9% confidence limits from these are 
shown together with the 0.1-1.0 keV power spectrum in 
Figure njtop panel). It is clear that both the 2193 and 
3510 sec periods are highly significant. To determine the 
errors on these periods we used the Cash statistic (Cash 
1979), which implies: 2193tgg and 3510tiQ4 sec - There are 
three other periods which are significant at > 99.9% level: 
1248ti5 sec , 1574±2? sec and 4844±^ sec. 

It is possible that some of these less significant periods 
could be produced as a result of aliasing between the most 
prominent periods (2193 sec and 3510 sec) and the window 
function or as a result of these two periods having a non- 
sinusoidal pulse shape. In order to investigate this further 
we have simulated a dataset with the same time points as 
the original data. However, this fake data only contains 
the two strongest signals (2193 sec and 3510 sec). The 
shape of these signals is obtained from the real data by 
using a simultaneous two-period fourier fit, where each of 
the two pulse shapes are represented with a second order 
fourier series. As a result we find that at least the 1574s 
period could be produced by the effects mentioned above, 
Figure (lower panel) . 



We discuss the possible origin of all of these periods in 
4. X-ray Spectral Analysis 

Taking the same time interval that the quasi-stable be- 
haviour is present as in the previous section, we phased 
events on the 2193 and 3510 sec quasi-stable periods re- 
ported in the previous section. Initially we extracted 6 
spectra covering different 2193 sec phase intervals. In ex- 
tracting these spectra we used apertures ~ 30 in radius 
centered on OY Car, chosen so that the aperture did not 
cover more than one CCD. This encompasses ~90 per- 
cent of the integrated PSF (Aschenbach et al 2000). A 
background spectrum was extracted from the same CCD 
on which the source was detected, scaled and subtracted 
from the source spectra. 

Since the response of the detectors is not well cal- 
ibrated at present below ~0.2keV, energies below this 
were ignored in the following analysis. The response file 
epn_ff20_sdY9_medium.rmf was used: this response in- 
cludes single and double pixel events. 

Ramsay et al (2001b) found that they could not well 
fit the integrated EPIC spectra of OY Car with a sin- 
gle temperature thin plasma model. A multi-temperature 
model gave a much better fit to the data. In their analy- 
sis Ramsay et al (2001b) concluded that a best fit for the 
integrated spectrum was obtained with a three tempera- 
ture plasma model with a partially covering absorber. The 
multi-temperature nature was supported from line H/He- 
like intensity ratios of Iron and Sulphur and also by the 
various line species seen in the high resolution RGS spec- 
trum. 

We fitted each of the 6 phase resolved spectra sepa- 
rately using the X-ray fitting package XSPEC. Each spec- 
trum was binned so that a minimum of 40 counts were 
present in each spectral bin. An initial analysis showed 
that a two temperature fit gave much better fits to the 
data. In these phase resolved spectra adding a third tem- 
perature component did not significantly improve the fit. 
Further, we find that the temperatures for each of ther- 
mal components were consistent in each of the 6 phase re- 
solved spectra (kT x =1.4keV and kT 2 =10keV). We there- 
fore fixed the temperatures of the thermal plasma mod- 
els at these temperatures. In each fit the free parameters 
were the absorption (we used a simple neutral absorption 
model) and the normalisation. The metal abundance was 
fixed at solar. 

The resulting fits are shown in tabulated form in Table 
^and in graphical form (along with the light curve folded 
on the 2193 sec period) in Figure^] The fits are reasonably 
good in the phase intervals 0.32-0.80, but in the remain- 
der it is less so (x^~1.5). In the case of the phase interval 
0.0-0.17 the fit was significantly improved when we fix the 
metal abundance at 1.5 solar (%J=1.35). It is clear from 
Figure 0] that there is no significant variability in the nor- 
malisations of the thermal plasma components. However, 
there is a variation in the column of the neutral Hydrogen 
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Fig. 4. From the top panel: using the time interval which 
showed evidence for a modulation on the 2193 sec pe- 
riod, we fold that data on that period in the 0.15-1.0keV 
energy band; the fitted neutral absorption column; the 
normalisation of the lower temperature component (fcTi); 
the normalisation of the higher temperature component 
(kT 2 ). 



density: it is anti-correlated with the folded intensity light 
in that it shows a maximum at intensity minimum. We 
show in Figure [5] the confidence contour plot in the (Nh, 
Tiow,norm) plane for the spectra covering the phase max- 
imum and minimum: there is a clear separation in the 
contours. 

We show in Figure El the X-ray spectra near 6-7keV 
for spectra taken from the phase maximum and minimum: 
the profile of the Fe Ka line does not differ significantly. 
This further strengthens the argument that the spectral 
changes over the 2193 sec period are most likely caused 
by the change in the amount of absorption. 

We also folded the data on the 3510 sec quasi-stable 
modulation: we find a similar result to that found for fold- 
ing on the 2193 sec. Namely, the temperature and normal- 
isation of the two thermal plasma components do not vary 
over phase, but the absorption component varies signifi- 
cantly. For brevity we do not show the results in Table or 
Figure form. 
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Fig. 5. Confidence contours in the (Nh, T\ , normalisation) 
plane for the phase minimum and maximum spectra. The 
contours correspond to 68%, 90% and 99% limits. 
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Fig. 6. The Fe Ka complex at the 2193 sec phase minima 
(0.5-0.6) (shown an asterix) and maxima (0.9-1.0) (shown 
as stars). 



5. Modelling the B band and UV data 

It is clear from Ramsay et al. (2001a) that the B and UV 
light curves exhibit clear modulation over the orbital pe- 
riod as well as asymmetric eclipse profiles. In this section 
we will describe and discuss the modelling of those data 
in detail. 

Our basic assumption is that all of the B and UV mod- 
ulation is caused by the changing visibility of different 
parts of (at least) partially optically and physically thick 
accretion disc and the eclipse of such disc by the secondary 
star. In our modelling we have taken the idea of 'fireflies" 
presented in Hakala, Cropper & Ramsay (2002) and de- 
veloped it further to suit the modelling of optically thick 
accretion discs. 
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Table 1. The results of phase resolved spectral fits folded 
on the 2193 sec period. The temperature of the cooler (Ti) 
plasma was fixed at 1.4keV, the hotter (T 2 ) at lOkeV. 



Hakala et al. (2002) use a swarm of emitting point 
sources (fireflies) to model the shape of an accretion 
stream in eclipsing AM Herculis systems. In their ap- 
proach the swarm is regularized to favour realistic 'ba- 
nana' shapes for the accretion stream (without fixing the 
exact location or shape of the stream). They then used a 
variant of a genetic algorithm to find the optimal accretion 
stream shape and location. 

If we want to apply a similar approach to locating the 
emitting regions in a disc accreting system, like OY Car, 
we need to make several modifications to the original code. 
In fact, the only assumptions that remain the same are the 
idea of 'fireflies' as a numerical model for the emission and 
the optimisation algorithm. We next discuss our model 
computation in detail. 

5.1. The model 

We first restrict the possible emission to originate within 
an accretion disc like volume around the primary. The 
boundary conditions for this volume are: maximum ra- 
dius 0.45a, maximum height ±0.02a. In addition we define 
a somewhat thicker disc like volume of maximum height 
±0.025a. This volume is assumed to have a homogenous 
optical thickness, that remains constant throughout the 
modelling. The reason for this is to provide estimates for 
the absorption from the different parts of the disc and thus 
guarantee an optically thick solution. 

To compute the optical light curve for any given dis- 
tribution of fireflies inside this boundary, the algorithm 
for the light curve generation procedes as follows. We first 
loop over all the orbital phases to be fitted and for each 
phase: 

— Select the fireflies that are not eclipsed by the sec- 
ondary; 

— Determine the absorption coefficient for each of the 
selected fireflies. This is done by computing the length 
of a path that light from any particular fly has to travel 
inside a predefined disc volume; 

— Sum up the emission from all visible flies at this phase 
(taking into account the variable absorption for differ- 
ent flies). 





0.85 0.90 



0.95 1.00 
Phase 



0.6 - 
0.4 - 

2 0.2 - 

o 

| 0.0 r 
o 

> -0.2 - 
-0.4 - 
-0.6 - 




X coordinate 



Fig. 7. Fits to the synthetic data. The whole light curve 
(top); zoomed in around the eclipse (middle); and the 
original (bottom, left) and reconstructed (bottom, right) 
brightness distributions for the fit. The bottom panel 
shows the ballistic stream trajectory and an accretion disc 
looking down on the orbital plane. 



The main difference between the 'firefly approach' and 
the 'normal approach' is that the firefly approach is (in 
analogue with hydrodynamics) a particle based solution 
whilst the 'normal' approach relies on fixed grid on which 
the solution is determined. The reasons behind our choice 
of the firefly approach here are the following ones. Firstly, 
one of the main axioms of the standard eclipse mapping 
routine (Home 1985) is that the eclipse profile of the disc 
is solely determined by the changing obscuration of the 
disc by the secondary. Thus outside the eclipse phases the 
light curve is constant. This implies that the same algo- 
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Fig. 8. Our best fit to the B-band data assuming no reg- 
ularization. The whole light curve is shown on top and we 
have zoomed in around the eclipse (middle). The result- 
ing firefly distribution together with a ballistic accretion 
stream trajectory is plotted in bottom (Primary at the 
origin, LI point at approx. coordinates of [0.7,0.0]). 



Fig. 9. Our best fit to the UV data assuming no regular- 
ization. The whole light curve is shown on top and we have 
zoomed in around the eclipse (middle). The resulting fire- 
fly distribution together with a ballistic accretion stream 
trajectory is plotted in bottom (Primary at the origin, LI 
point at approx. coordinates of [0.7,0.0]). 



rithm could not be used to fit the rest of the light curve. 
Normally the out-of-eclipse variation, often seen in quies- 
cent state CV light curves, has to be removed before the 
fitting, which can cause undesirable bias in the results. 
Rutten (1998) developed the eclipse mapping by includ- 
ing a 3D shape for the disc. His model could be used for 
out-of-eclipse light curve modelling, but it has it's limita- 
tions as well, i.e. we have to know and fix the 3D shape of 
the disc beforehand. 

The advantage of the firefly modelling is the lack of 
any grids. The flies are free to move within a predefined 



volume and their final distribution represents the location 
or distribution of emission within the disc. However, our 
approach does require that the disc has an axisymmetric 
homogenous density structure, as this is assumed for the 
computation of optical thickness effects that dominate the 
out-of-eclipse light curve shape. 

5.2. Model Regularization 

As in most of inversion problems, there are an infinite 
number of solutions if no regularization is employed. In 
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Fig. 10. Our best fit to the B-band data assuming axisym- Fig. 11. Our best fit to the UV data (with axisymmetric 

metric regularization. The whole light curve is shown on regularization). The whole light curve is shown on top 

top and we have zoomed in around the eclipse (middle), and we have zoomed in around the eclipse (middle). The 

The resulting firefly distribution together with a ballistic resulting firefly distribution (bottom), 
accretion stream trajectory is plotted in bottom (Primary 
at the origin, LI point at approx. coordinates of [0.7,0.0]). 



order to overcome this, we have (like many previous au- 
thors) decided to search for the most axisymmetric image 
that can reproduce the observations (axisymmetric about 
the accretion disc rotation axis). In the case of the firefly 
approach this simply means a regularization term, where 
we have defined the disc in i sectors, and counted the num- 
ber of flies in each sector Ni (here we have chosen to use 36 
sectors). Then we maximize the entropy of such division 
i.e. we first compute 



Pi = 



(i) 



where pjis the fractional flux in each sector and the en- 
tropy term is 



36 



S = 1o s(k) 

i=l 

finally our merit function to be minimized is 



(2) 
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Fig. 12. The radial distribution of fireflies in B (solid) and 
UV bands (dashed). The hot spot in B band is visible at 
a distance of 0.3a. 

F = xl - \S (3) 

where xt is the standard reduced x 2 term, A is the 
Lagrangian multiplier and S the entropy term. 

5.3. The fitting algorithm 

We carry out the actual fitting using the same genetic 
algorithm variant as in Hakala et al. (2002). The only 
difference is that in this case the initial population of so- 
lutions for fitting consists of disc like swarms of random 
fireflies. These are then evolved in similar fashion to fit the 
eclipse profile and the overall light curve simultaneously. 
As the whole light curve has 532 points in B and 514 points 
in UV, we have decided to emphasize the eclipse phases 
(0.95-1.05) by summing the Chi-square values from this 
phase range three times. We typically use 50 individual 
swarms in each generation and after around 1500 genera- 
tions, the fits stop evolving in a reasonable time. 

In order to test our approach we have created a syn- 
thetic light curve using a known disc brightness distri- 
bution, added appropriate levels of noise and fitted that. 
The results of this exercise are shown in Figure The top 
panel shows the whole orbital light curve together with the 
fit, the middle panel the phase range close to the eclipse 
phase and the bottom panel the synthetic brightness dis- 
tribution (left) and the reconstructed distribution (right). 
The fits to the synthetic data are encouraging, with the 
reconstructed disc being close to the input disc. There is, 
however, a small amount of smearing, which is most likely 
caused by finite phase resolution, regularization and/or 
noise in the data. 

5.4. The Results 

We have plotted in Figures|Hland|nithe resulting brightness 
distribution from a fit to the B and UV band OM data 



in a case where no regularization was applied. The fit is 
excellent, but the resulting brightness distribution suffers 
slightly from the 'scissors' effects, i.e. the point spread 
function of the map has a scissors shape. This is because 
there is not enough information in the data to uniquely 
locate point like emission. Thus point source emission is 
spread out along the lines that mark the shadow of the 
secondary at the ingress and egress phases. This is a well 
known effect and was demonstrated by Home (1985). 

Using the regularization described earlier we obtain 
the fits for B and UV bands shown in Figures ED and 
^^respectively. The coordinate system for the fireflies is a 
Cartesian one, where the primary star is located in the ori- 
gin and the X-Y plane is the orbital plane. The coordinates 
of the LI point are roughly (0.7,0,0) and a ballistic freefall 
stream is overplotted. The observer lies 'above' the orbital 
plane, i.e. she/he would have a positive Z-coordinate. 

Figure IT2"1 shows the radial histograms of the distribu- 
tion of flies. The resulting brightness distributions seem to 
indicate that whilst the UV emission comes mostly from 
the white dwarf or the boundary layer and the inner disc, 
the optical B band emission also contains a faint com- 
ponent from the accretion stream - disc interactions, i.e. 
the hot spot (and maybe even from the freefall part of the 
stream). This is also rather obvious directly from the light 
curve shapes, given the difference in eclipse assymmetry 
between the two passbands. 

Furthermore, we find that the hot spot is located ex- 
actly on the ballistic trajectory, as expected. The distance 
is ~ 0.3a from the white dwarf, which confirms the val- 
ues reported by Wood et al. (1989). In addition to this 
we find marginal evidence for optical emission from the 
accretion stream prior to the main hot spot impact re- 
gion. Our results show, that in the B band 90% of the 
emission comes within a distance of 0.29a from the white 
dwarf and 95% within 0.33a. This implies that the hot 
spot contributes roughly 5% of the total B band emission. 
The same radii in UV band are slightly greater in compar- 
ison (0.32a and 0.38a). At face value this implies that the 
B band radial emission distribution is concentrated more 
towards the white dwarf than in UV passband. This differ- 
ence could be due to the fact that whilst our UV data was 
taken 4 days after an outburst, the B band data was taken 
more than a month later. Dwarf novae discs are thought to 
be more uniformly bright during the outbursts. Although 
the V band mag was similar in both observations this does 
not preclude a difference in the size of the accretion disc 
in the two observations. 

6. Discussion 

6.1. The origin of the X-ray modulation 

Ramsay et al (2001a) found that in X-rays the duration of 
the eclipse ingress and egress is 20-30sec, which is roughly 
the same time as that expected for the secondary star to 
eclipse the diameter of the white dwarf. These authors 
concluded that the X-ray emission region is hence ex- 
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tended in white dwarf longitude, probably centered near 
the spin equator. However, Wheatley & West (2003), us- 
ing a re-analysis of the same data, conclude that the X-ray 
emission originated close to the upper pole on the white 
dwarf. In the latter scenario, an accretion column, or cur- 
tain, of material will rotate on the spin period of the white 
dwarf. 

This is consistent by our spectral modelling in 21which 
indicates that the most likely cause of these quasi-stable 
modulations is a variation in the absorption column over 
this period. If the 2193 sec period was the spin period 
of the white dwarf, how would that manifest itself as a 
variation in the absorption column? In the case of inter- 
mediate polars, some systems show an energy dependent 
X-ray variation at the spin period. This has been seen in 
EX Hya (Rosen, Mason & Cordova 1988, Cropper et al 
2002) and is thought to be due to observing the emis- 
sion site through an accretion curtain whose absorption 
column varies as the white dwarf rotates. 

A second possibility may be due to a variation in the 
distance from the inner edge of the accretion disc to the 
white dwarf. Further, this distance could vary in azimuth 
in the binary frame. Thus the white dwarf-accretion disc 
interaction would lead to a variation at a beat period and 
could account for why the modulation is only seen for 
around 1/3 of the observation duration. This hypothesis 
requires that the 2193 sec period is the spin period of the 
white dwarf (u>), whilst the 3510 sec modulation is the beat 
period between the spin and the 5436 sec orbital period 
(fi). The beat period based on the 2193 s spin period 
would be (u> — ft) 3676 sec, which is 6.6er distant from the 
best fit period. However, if the spin period is actually 1.7a 
shorter than the best fit value (giving 2133 sec) then the 
observed 3510 sec period would be exactly the expected 
lo — Q, beat period. In this scenario, any very short periods 
(for instance the 17.94 and 18.15 sec periods reported by 
Marsh & Horne, 1998) would have to be explained in other 
terms, for instance white dwarf oscillations. 

Finally, Marsh & Horne (1998) suggested that one of 
the ~18 sec periods could represent lo, or alternatively 
that one of the periods was the Keplerian period and the 
other was the beat period between this and the spin pe- 
riod (implying a spin period of ~1550 sec). If ~1550 sec 
was the spin period, then most of the amplitude peaks 
in our Figure ^ could be explained as naturally occuring 
sidebands of the spin and orbital frequencies in the fol- 
lowing manner: 1248 sec (lo + fl), 1574 sec (to), 2193 sec 
(lo - £1), 4844 sec (2(u - ft)). The main drawback with 
this hypothesis is that there is no simple explanation for 
the very significant 3510 sec period. Also, it is possible 
that not all of the 5 periods mentioned above are in fact 
physical (cf §3). 

Recently, the short period variability in cataclysmic 
variables has been extensively discussed in a series of 
papers by Warner & Woudt (most recently in Warner, 
Woudt & Pretorius (2003)). In conclusion they have found 
that typically Pqpo ~ 15-Pcwo- In addition they also find 
second type of DNO's with: Pi p dno ~ ^Pdno- Assuming 



that one of the periods seen here, 2193s or 3500s, would 
be a QPO period the above mentioned relationship would 
then predict DNO periods around 150s and 230s respec- 
tively. There is no evidence for such periods in our data. 
However, this is expected since the DNO and QPO peri- 
ods are normally found during outbursts, and during our 
observation the source was in quiescence. 

We conclude that based on our data, we cannot provide 
a unique explanation for the 2193 sec and 3510 sec periods. 
We can, however, conclude that these periods occur due 
to changes in the line of sight absorbing column. 

6.2. The Optical Light Curve Modelling 

Much previous work has been done on the eclipse pro- 
files of OY Car ranging from infrared to X-rays (see for 
instance Berriman 1984, Cook 1985, Krzeminski & Vogt 
1985, Schocmbs 1986, Berriman 1987, Wood et al. 1989, 
Rutten et al. 1992, Horne et al 1994, Bruch, Beele and 
Baptista 1996, Pratt et al. 1999, Ramsay et al. 2001a, 
Wheatley & West 2003). Wood et al. (1989) use the opti- 
cal eclipse profiles to measure accurate system parameters 
and assuming a ballistic free fall trajectory for the accre- 
tion stream they deduce a disc outer radius of approxi- 
mately 0.31a. Rutten et al. (1992) re-analysed the set of 
eclipse profiles from Vogt (1983) that cover the rise into 
an outburst. Their eclipse mapping analysis shows that 
whilst the quiescence eclipse profiles are dominated by the 
white dwarf and the hot spot, the rest of the disc becomes 
dominant in the outburst. The HST UV observations by 
Horne et al. (1994) demonstrate that the eclipse profiles 
become more symmetric and 'tophat shaped' towards the 
shorter wavelengths. This indicates that the white dwarf 
and/or the boundary layer is the major source of emission 
in the UV, whilst the hot spot contribution rises towards 
the optical wavelengths. 

Our light curve modelling confirms the previous results 
on the origin of the quiescent eclipse profiles. In particular 
they indicate, that both in B and UV bands, the majority 
of the emission arises close to the white dwarf, possibly 
either the boundary layer or the white dwarf itself. This 
is also supported by the ingress and egress time analysis 
of the the same B, UV dataset by Ramsay et al. (2001a). 

The freefall accretion stream and especially the hot 
spot, where the stream hits the disc outer edge is mainly 
emitting in the B band and hardly visible in the UV maps. 
Our modelling also suggests that in addition to the hot 
spot, the ballistic accretion stream might contribute to 
the B band emission as well. 

Our novel approach allows simultaneous fitting of the 
eclipse profiles and orbital light curves of cataclysmic vari- 
ables. This removes bias introduced by artificially having 
to remove the orbital modulation prior to the eclipse pro- 
file mapping. Furthermore, the 'firefly approach' has po- 
tential to be easily converted into full 3D mapping, if a 
proper regularisation can be employed. 
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7. Conclusions 

We have reanalysed the XMM-Newton observations of the 
dwarf novae OY Car. In particular, we have reanalysed 
the X-ray light curve and spectrum and developed a new 
model to fit the OM B band and UV data. 

We have investigated the X-ray light curve over the 
time interval where the previously reported 2240 sec mod- 
ulation was strongly seen. This has allowed us to refine 
this period to 2193 sec. We also detect a modulation at a 
period of 3510 sec which we also believe to be significant. 
We have performed phase resolved X-ray spectroscopy on 
these periods and find that the soft X-ray light curve is 
anti-correlated with the absorption column density. This 
implies that the light curve modulation is due to us observ- 
ing the X-ray emission region through a rotating accretion 
column or curtain. Although, not conclusive, this supports 
the hypothesis of Ramsay et al (2001a) that the 2193 sec 
is the spin period of the white dwarf. Furthermore, this is 
in agreement with findings of Wheatley & West (2003), 
where they conclude that the X-ray emission probably 
originates near the polar region on the white dwarf sur- 
face. 

To model the OM B and UV light curves, we have 
applied the fireflies model of Hakala et al (2002). We sim- 
ulated light curves using a given accretion disc brightness 
distribution and used our technique to fit the synthetic 
light curves. The resulting distribution of fireflies is en- 
couraging. We then applied the model to the OM data and 
find that that the hot-spot, where the accretion stream 
hits the disc, is exactly on the ballistic trajectory. Further, 
we find some evidence that the B band emission is more 
concentrated than the UV emission, although this maybe 
due to the fact that the B band data were taken 1 month 
after the UV data. Our technique can be applied to other 
disc accreting systems. 
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